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Abstract 

The Brownian loop soup introduced in |5j is a Poissonian realization from a cj-finite 
measure on unrooted loops. This measure satisfies both conformal invariance and a 
restriction property. In this paper, we define a random walk loop soup and show that 
it converges to the Brownian loop soup. In fact, we give a strong approximation result 
making use of the strong approximation result of Komlos, Major, and Tusnady. To 
make the paper self-contained, we include a proof of the approximation result that we 
need. 



1 Introduction 

The Brownian loop soup with intensity A, which we define below, is a Poissonian realization 
from a particular measure on unrooted loops in C that satisfies both conformal invariance 
and a property called the restriction property. A realization of the loop soup consists of a 
countable collection of loops. In a fixed bounded domain there are an infinite number of 
loops that stay in D\, however, the number of loops of diameter at least e in the bounded 
domain is finite. A corollary of conformal invariance is scale invariance: if ^ is a realization 
of the Brownian loop soup and each loop is scaled in space by 1/iV and in time by l/N"^, 
the resulting configuration also has the distribution of the Brownian loop soup. In this 
paper, we will show that the Brownian soup is a limit of random walk soups. There are 
two natural approaches to showing this. One is a "weak" limit to show that the Brownian 
loop measure is a weak limit of random walk measures (this requires some care since the 
measures are infinite). However, we choose the more direct "coupling" approach of defining 
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the random walk loop soup and the Brownian loop soup on the same probability space so 
that the realizations are close. Since a realization is a countable collection of loops, it is a 
little tricky to say what it means for the realizations to be close. We will prove, in fact, 
that in a bounded domain except for an event of small probability, there is a one-to-one 
correspondence between the Brownian loops and the random walk loops if we restrict to 
loops that are not "too small" . The Brownian loops and random walk loops that correspond 
with each other will be very close. We will use the dyadic approximation scheme as in [2] to 
establish the strong approximation of the two soups. 

We start by defining the Brownian loop measure. It is easier to define the loop measure 
first on rooted loops. A (rooted) loop is a continuous function 7 : [0,t-y] C with 7(0) = 
7(t^). We will only consider loops with < < 00. Let C denote the set of all loops and Ct 
the set of loops 7 with t^ = t and 7(0) = 7(^7) = 0. The Brownian bridge measure /x'^'' is the 
probability measure on loops induced by a Brownian bridge , i.e., by := Wt — tWi, < t < 
1, where Wt is a standard two-dimensional Brownian motion. The measure is supported 
on Ci. The (rooted) Brownian loop measure is the measure on Ci x C x (0, 00) given by 

/X = fi^'^ X area x ^ dt]- 

The measure /i induces a measure on C, which we also denote by fj,, by the function (7, z, t) 1-^ 
7, where 7 is 7 scaled (using Brownian scaling) to have time duration t and translated to 
have root z. In other words, 

7(s) = z + t^/^7(s/t), 0<s<l 

This measure is clearly translation invariant, and it is straightforward to check that if r > 0, 
then n is invariant under the Brownian scaling map (7,2;, t) 1-^ {'y,rz,r'^t). 

Let us denote by fi^^{z) the probability measure on loops induced by a Brownian bridge 
of time duration t rooted at z. Then the measure /i (as a measure on C) can be written as 

An unrooted loop is an equivalence class of (rooted) loops under the equivalence 7 ~ Or'j 
for every r G R, where 6r'~f{s) = 7(3 + r) (here we consider a rooted loop 7 of time duration 
t^ as a continuous function 7 : M ^ C with 7(5 + 1^) = 7(5) for all s). The unrooted loop 
measure /i" is the measure obtained from jj, by "forgetting the root." A rooted Brownian 
loop soup with intensity A is a Poissonian realization from A/x. An (unrooted) Brownian loop 
soup is a realization from A/x". One can obtain an unrooted loop soup by starting with a 
rooted loop soup and forgetting the root. 

If D is a. domain in C we let fi£),fJ'^D denote /x, /i" restricted to loops that lie in D. 
The family of measure {/x"/)} clearly satisfy the restriction property, i.e., ii D' G D then 
/x"^,/ is /x"^ restricted to curves lying in D'. It is also shown in 3] that the family satisfies 
a conformal invariance property, i.e., ii f : D ^ D' is a. conformal transformation, then 
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/ ° /^"d = A^^D' 5 if the quantites are suitably interpreted. In particular, if 7 is a curve lying 
in we define / o 7 to be the curve in D' , reparametrized by the conformal map; see jS] for 
details. The measure /i on rooted loops is not conformally invariant. 

In this paper we study the loop measure for simple (nearest neighbor) random walks on 
the integer lattice Z^, which we can consider as a subset of C. The rooted loop measure ^™ 
gives each (nearest neighbor) random walk loop in 1? of length 2n measure {2n)^^ 4~^". The 
unrooted loop measure is obtained from the rooted loop measure by "forgetting the root" . 
It is almost the same measure as that obtained by giving measure 4"^" to every unrooted 
loop of length 2n. (If a loop of length 2n is obtained by taking a loop of length n and 
repeating the same loop again, then this unrooted loop does not get full measure 4~^" under 
our random walk loop measure; these exceptional loops are an exponentially small subset 
of the set of all loops so it is not important whether we give these unrooted loops measure 
^-2n (1/2) 4~^".) We will focus on the rooted measure in this paper. A rooted random 
walk loop of length 2n can also be considered as a continuous path 7 : [0, 2n\ — >■ C by linear 
interpolation. We will call a Poissonian realization from A /x™ a rooted random walk loop 
soup (with intensity A). 

In this paper, we make a precise statement that the random walk loop soup, appropriately 
scaled, approaches the Brownian loop soup. We will define {A\, A\) on the same probability 
space so that is a realization of the Brownian loop soup with intensity A, and is a 
realization of the random walk loop soup with intensity A. We consider the loops in ^a as 
curves in C by linear interpolation. Note that ^a is a (random) countable set of curves and 
^A is a (random) multi-set (i.e., a set where some elements can appear more than once) 
of lattice curves. For each positive integer N we define ^a.a^ to be the collection of loops 
obtained from ^a by scaling space by More precisely, 

Ax,N = {$iv7 : 7 e Ax}, 

where = t^/N"^ and 

$^7(t) = Ar-i7(tAr2), Q<t<tjN^. 

Note that the scaling rule implies that ^A,Af is a realization of the Brownian loop soup with 
parameter A. We define 

Ax,N = {^Ni ■■ 7 e Ax}, 

where t^^^ = t^/{2N^) and 

|.^7(t) = N-^-f{t2N^), 0<t< t^/{2N^). 

The scaling is slightly different for the random walk loops because the covariance of a simple 
two-dimensional random walk in 2n steps is nl as opposed to 2nl for a Brownian motion 
at time 2n\ roughly speaking, this is because in 2n steps, the random walk moves about n 
steps horizontally and n steps vertically. 
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We will prove the theorem below. The ideas in the proof are simple and flexible. However, 
due to discretization, stating the result is a little unwieldy. To aid, we introduce the following 
auxiliary functions. For t > (5/8)A^~^ and positive integer k, we let 

Also, for z E C, zq E Zi"^ we define 

Mz) = ^ if max{|Re{Nz-zo}|,|Im{Nz-zo}|}<i. 

The definition of ip^ if Nz happens to fall on a bond of the dual lattice of is irrelevant 
for our theorem. 

Theorem 1.1. One can define on the same probability space Ax and Ax such that: 

• For each A > 0, is a realization of the Brownian loop soup; the realizations are 
increasing in A. 

• For each X > 0, Ax is a realization of the random walk loop soup; the realizations are 
increasing in X, 

and such that the following holds. Let Ax,n O'nd Ax,n be as defined above. Then there exists 
a c > such that for every r > 1,N,X and every 2/3 < 9 < 2, except perhaps on an 
event of probability at most c(A + l)r'^ N"^'^^ , there is a one-to-one correspondence between 
{7 eAx,N : h > N<>-\ |7(0)| < r} and {7 G Ax^n : ^fNit,) > N<>-\ |^^(7(0))| < r}. // 
7 G Ax,N and 7 G Ax,N are paired in this correspondence, then 

\t^-t^\ < 5/8N~^ 

sup |7(st^) -7(s^)| < cA^^MogA^. 

0<s<l 

The outline of the paper is as follows. In the next three sections, we define the random 
walk loop soup, state the strong approximation result between random walk bridges and 
Brownian bridges that we need, construct the probability space on which both the random 
walk and Brownian loop soups are defined, and verify that the construction satisifes Theorem 
11.11 The next section concerns the soups in bounded domains. Here we establish a similar 
result to the theorem above, although the error terms are somewhat larger. The remainder 
of the paper gives a self-contained proof of the strong approximation result that we need. 

2 Random walk loops and random walk soup 

A (rooted) random walk loop of length 2n in (which will be considered as a subset of C) 
is a (2n + l)-tuple u = [uq, . . . ,U2n] with \ujj — ujj^i\ = 1 and uq = U2n- A loop can be 
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considered as a curve 7 : [0, 2n] C; here 7(m) = uJm for integer m and 7(t) is defined for 
other t by hnear interpolation. Let denote the set of random walk loops of length 2n 
and be the set of such loops with uq = z. There is a natural one-to-one correspondence 
between £° and given hj u ^ z + lj. We let C = U„>i£„, = U„>i£^. We let u 
denote the random walk loop measure C, i.e., the measure that assigns measure 4~^" to each 
u G C2n- We write for u restricted to C^. For fixed z, it is straightfoward to show 
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The second equality follows from the fact that the probability that a two-dimensional simple 
random walk returns to the origin at time 2n is the square of the probability that a one- 
dimensional simple random walk returns to the origin; see (j2)) below to see why this is true. 
The final equality is derived from Stirling's forumla with error: 
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We define the (rooted) random walk loop measure fi™ to be the measure that assigns 
measure {2n)~^i'{uj) = (2n)~^4~^" to each uj G Cn- A rooted random walk loop soup of 
intensity A is a Poissonian realization from the measure A /i™. We will obtain a realization 
by first using a Poisson point process to select a multi-set of ordered pairs {n, z), where the 
length of the loop is 2n and it is rooted at z. Then, given (n, z), we choose a loop from the 
appropriate random walk bridge measure. In other words, we write 
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2n 
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Let 



N{n,z]t), n G {1,2,3,...}, z G Z^, 
be independent Poisson processes (in the variable t) with parameter 



2n ^ 
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Let 



L(n, z\rn), n E {1, 2, 3, . . .}, z G Z^, m G {1, 2, 3, . . .} 

be independent random variables, independent of the N{n,z]t), taking values in C^; the 
distribution of L{n, z; m) is the probability measure of a random walk bridge of 2n steps, 
i.e., which is the uniform probability measure on If G let 

N{n,z;t) 

Jt{^) = ^ l{L(n,2;fc)+2=a;}- 
fc=l 
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Note that {Jt{uj) : u; G £} is a collection of independent Poisson processes; the process Jt{uj) 
has parameter (2r?,)^^4^^" if G £„. We could equally well have constructed the loop soup 
starting with these Poisson processes. We have chosen the longer construction because it 
will be useful for coupling the loop soup with the Brownian loop soup. 

Although we have used t for the time parameter of the Poisson processes, by choosing 
t = A we get an increasing family of realizations of the loop soup A\ parametrized by A. We 
think of the loop soup of intensity A as a multi-set Ax of loops where loop u appears Jx{uj) 
times in A\. 



3 Strong Approximation 

If Sn denotes a simple random walk, we define St, < t < oo, by linear interpolation. 
The key to the coupling is the following result, due to Komlos, Major, and Tusnady, which 
shows that a simple random walk bridge and a Brownian bridge can be coupled very closely. 
Because the form of the result we need is slightly different than that proved in [2], we have 
included a proof in the final section. By a simple random walk bridge of time duration 2n 
we will mean a process Xt, < t < 2n, that has the law of St, < t < 2n conditioned to 
have S2n = 0. 

Lemma 3.1 (Dyadic approximation). There exists a c < oo such that for every positive 
integer n, there exists a probability space {Q, T , P) on which are defined a one- dimensional 
Brownian bridge Bt, < t < 1 and a one- dimensional simple random walk bridge Xt, < 
t < 2n such that 



P{ sup \{2n)-^^^X2ns - Bs\ > crT^I'^ logn} < cn 



30 



0<s<l 

Proof. This is a special case of Theorem 16.41 □ 

Remark. The choice of 30 as the exponent on the right-hand side is arbitary. The same 
result holds with error 0{n~^) for any r > 0, with suitably chosen c = c^. 

Corollary 3.2. There exists a c < oo such that for every positive integer n, there exists a 
probability space P) on which are defined a two-dimensional Brownian bridge Bt,0 < 

t < 1 and a two-dimensional simple random walk bridge Xt, < t < 2n such that for each 
n, 

P{ sup \n-^/^X2ns -Bs\> cn-^l'^ logn} < cn''^^ . 

0<s<l 

Proof. If Sj, Sj are independent one- dimensional simple random walks, then 

._ S} + iS] 

is a two-dimensional simple random walk (written in complex form). Conditioning on S2n = 
is the same as conditioning on S'g^ = = 0. In other words, we can obtain a two- 
dimensional random walk bridge as the product of two independent one-dimensional random 
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walk bridges. Hence we can construct the probability space as the product of two probability 
spaces as in the lemma. □ 

The following corollary is in the form that we will need in the rest of the paper. We have 
chosen to write it out in detail because of this. Recall that if a process is defined only for 
integer times, we extend its definition to non-integer times by linear interpolation. 

Corollary 3.3. There exists a c and a probability space {Q, JF, P) on which are defined process 
Qn,z,m^ gn,z,m^ n = 1, 2, . . . , z G Z^, m = 1, 2, . . . such that 

• the processes 

^n,z,m^ 0<t<l; n = 0,l,...; zEZ^, m = l,2, 
are independent two-dimensional Brownian bridges; 

• the processes 

gn,z,m^ r;, = 0, 1,...; z E , j = 0, . . . ,2n, m = l,2, 

are independent and Sj'^'"^,j = 0,...,2n has the distribution of a two-dimensional 
simple random walk conditioned so that 5*2^^'™' = 0. 

P{ sup |n-^/2^2Y;" - ^."'''"1 > cn-'/^ logn} < cn-^'' . 

0<s<l 

Proof. Take products of the probability spaces in the previous corollary. □ 

For the remainder of this paper we fix the probability space (fi, JF, P) as in the previous 
corollary. On this probability space are defined the independent valued random variables 
L(n, z] m) as in Section |21 

f (rn y. rn\ — qn,z,m qn,z,m-\ 

Also, we have independent, identically distributed random variables L(n, z; m) taking values 
in Ci. 

L(n, z- m) = S"'"'"^ 

These have the distribution of a two-dimensional Brownian bridge of time duration 1. We 
assume that these are coupled as in the proposition. 
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4 Constructing the Brownian loop soup 



In this section we will show how to construct the Brownian loop soup in a way that is highly 
correlated with the random walk loop soup. We will restrict to the rooted Brownian loop 
soup restricted to loops of time duration at least 5/8; to get a complete realization one can 
attach an independent realization of the loop soup with loops of time duration less than 
5/8. These small loops will not be coupled with the random walk loops. Let N{n, z; t) be a 
collection of independent Poisson processes with parameter 

"+(5/8) _ds_ ^ 1 ^ ^ l_ 

„_(3/8) 27rs2 ~ 27r(n + (5/8)) (n - (3/8))) ~ 2Txn'^ Snn^ ^ n^'' 

Recall the definition of N{n,z;t) from Sectional and note that qn — Qn = O(^) . (We 
have chosen to couple random walk loops of 2n steps with Brownian loops of time duration 
ra — (3/8) to n + (5/8). The particular choice of interval [—3/8, 5/8] was chosen so that g„ and 
g„, agree up to an error of size 0{n~^).) It is easy to see that we can couple N{n, z; t), N{n, z, t) 
on the same probability space so that: 

• {N{n, z;t)} are independent Poisson processes with parameter 

• {N{n, z;t)} are independent Poisson processes with parameter g„, ; 

• There is a c such that for all n, z, t, F{N{n, z; t) ^ N{n, z;t)} < t |g„ — < ctn~'^. 

In fact, we can let N{n, z;t),n = l,2,...,z G Z^, be independent Poisson processes with 
parameter 1 and then set 

N{n, z; t) = N{n, z; qj), N{n, z; t) = N{n, z, ; qj). 

Assume without loss of generality that on this probability space we have independent copies 
of the coupled processes 

L{n,z;m), L{n,z;m) 

as in Sectional independent complex- valued random variables Y{n, z; m) that are uniformly 
distributed on the square {x+iy : |x| < 1/2, \y\ < 1/2}, and independent real-valued random 
variables T{n, z; m) with density 

^' n < s <n ^ — . 



{n + l){n - I) 3 ,5 

^2 



We construct the rooted Brownian loop soup (restricted to loops of time duration at least 
5/8) as follows: 

• N{n, z\ t) will be the number of rooted loops that have appeared by time t whose root 
is in the unit square centered at z and whose time duration is between n — (3/8) and 
n-F(5/8); 
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• scale the bridge (of time duration 1 and rooted at 0) L{n,z,m) so that it has time 
duration T{n, z, m); and then translate it so that its root is z + Y{n, m); we call this 
final loop L*{n, z,m). 

Then it is easy to see from the definition that the collection of loops 

Ax = {L*{n, z] m) : N{n, z; X) > m, n e Z+, z G Z^} 

is a realization of the Brownian loop soup with intensity A (restricted to loops of time 
duration at least 5/8). We can then extend Ax to a realization of the Brownian loop soup 
by adding an independent realization of loops of time duration less than 5/8. Recall from 
the discussion in the Introduction that ^a.a^ is also a realization of the Brownian loop soup. 
On this space we also have the scaled random walk soup Ax,n- 

We will now show that this coupling satisfies the conclusions of Theorem 11.11 Without 
loss of generality, we may assume that Ar^iV^"^^ < 1; in particular, A < A^^. First, note that 

F{N{n, z; A) ^ N{n, z; A) for some <n <oo, \z\ < rN} 

\z\<rNn>Nf> 

Hence, except for an event of probability 0{Xr'^N'^~^^), N{n,z]X) = N{n,z]X), for n > 
N^, \z\ < rN. Also, 

¥{N{n,z;X) > for some iV^/^ <n<N^, \z\ < r N} < cr^ N^¥{Y > N^} < cr^N'^, 

where F is a Poisson random variable with expectation cN^. The last estimate uses an easy 
estimate on Poisson random variables; in fact for positive integers N, F{Y > N^} < F{Y > 

Let 

Then Z is Poisson with 

\z\<rNn>N^ 

Hence, F{Z ^ 0} < cAr^A^"^ < cAr^A^^-se^ Therefore, 

F{N{n,z,X) > for some n > A^^ \z\ < rN} < cXr^N'^~^^ 

and similarly the same estimate holds with N{n,z,X) replacing N{n,z,X). 

Let us denote the loops L*{n, z; m) and I/(n, z;m) + z by 7^,2,™ and 7n,2,m- Let A = A^^r 
be the event 

A = { sup \'yn,z,mist^) " ln,z,rn{st^)\ > Ca log A^*' 
0<s<l 
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for some \z\ < rN, N'^'^ < n < N^,m < N^}. 
Here we use C2 for the constant c from Corollary 13.31 Then the corollary tells us that 

On the intersection of A'^ and the events that 

N{n, z; A) = N{n, z;X), n> N\ \z\ < rN, 

N{n, z; A) < N^^^ <n<N^, \z\ < rN, 

N{n, z; A) = N{n, z,\) = 0, n> N^, \z\ < rN, 
the coupling satisfies the conclusion of Theorem 11.11 

5 Bounded domains 

In this section, we let D denote a simply connected domain in C containing the origin, 
contained in the unit disk B = {z e C : \z\ < 1}. li e > 0, let = {z e D : dist(2;, dD) > e}. 
Recall the definition of the loop measure /i from the Introduction. We say that a loop 7 is 
in D if 7[0,t^] C D. 

Proposition 5.1. There is a c < 00 such that «/0 < e < tg , then the n measure of the set 
of loops in D of time duration at least to that are not in is bounded above by cetQ 

Proof. The measure we are interested in is given by 

Li i^'i'-)""^. 

where pt{z) = P{0 < dist(-B^'' [0,t],(9D) < e} and B^"^ denotes a Brownian bridge of time 
duration t rooted at z. By time reversal, we can see that, 

P{0 < dist(Sf [0,t],aD) < e} < 2P{0 < dist(S^'^ [0, -],aD) < e,B^' [0,t] C D}. 

But, 

P{0 < dist(S^'' [0, ^], dB) < e, B"^' [0, t] C ©} = 

2t t 
lim -r^rP^iO < dist(5[0, -1, dB) < e, B\0, t] C D, \Bt - z\ < S], 

where Bf denotes a standard Brownian motion. By the strong Markov property and the 
"gambler's ruin" estimate^, 

P"{0 < dist(5[0, 83) < e; B[0, h] C3}< cet^^/^ 

^The gambler's ruin estimate states that the probabiUty that a one dimensional standard Brownian motion 
starting at e > stays positive up to time t is bounded above by cet~^^^. 
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Given this event, the probabihty that \Bt — z\ < 6 is bounded above by c6'^/t. Hence 

Mz)<ce6t-'/\ 

and the result follows by integrating. □ 
Proposition 5.2. Suppose D is a simply connected domain contained in the unit disk. There 

5 /4 

is a c < oo such that if < e < t^ , then the /i measure of the set of loops in D of time 
duration at least to that are not in is hounded above by c^l"^ to'^^ . 

Proof. The measure we are interested in is 

f oo 

Pt{z) dt dz, 




D J to 



27rt2 



where pt{z) = P{0 < dist(-B^'' [0,t],dD) < e} and B^"^ denotes a Brownian bridge of time 
duration t rooted at z. By time reversal, we can see that. 



P{0 < dist(5^^'' [0, t], dD) < e} < 2P{0 < dist(5,^'^ [0, dD) < e, B^^' [0, t] C D}. 



But, 



< dist(5,^^ [0, dD) < e, B^/ [0, t]cD} = 

2t t 
lim — Pno < dist(5[0, -1, dD) < e, B\0, t] C D, \Bt - z\ < 6}, 
s->o+ 0^ ' ' 2 

where Bt denotes a standard Brownian motion. By the strong Markov property and the 
Beurling estimate (see Lemma 1^31 below) . 

P"{0 < dist{B[0,^],dD) < e;B[0,h] C D} < ce^^H-^/^ 

Given this event, the probability that \Bt — z\ < 6 is bounded above by c6^/t. Hence 

p,{z)<ce'/H-'/\ 

and the result follows by integrating. □ 

Lemma 5.3 (Beurlng estimate). Let Bf denote a standard two-dimensional Brownian 
motion. There is a c such that if'-/ : [0, oo) ^ C is any continuous curve with \'j{0)\ = r and 
lim^^oo |7(^)| = OO; then 

F'{B[0, t] n 7[0, cx)) = 0} < c (^-^^ \ \z\< r. 
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Proof. If t is replaced by at = inf{s : |-Bs| = then this lemma follows from the Beurling 
estimate which is a corollary of the Beurling Projection Theorem Theorem V.4.1]. We 
will assume that estimate and show how the estimate for fixed times t can be deduced from 
the result for the stopping times at- By scaling it is enough to do t = 1. Let Tn = inf{t : 
\Bt-Bo\ = 2-"}. We will show that F{B[0, 1 A Tq] n 7[0, oo) = 0} < cr 5 . 
First, note that 

P"{S[0,1 Aro] n7[0,oo) = 0} 

< ¥'{B[0, To] n 7[0, oo) = 0} + W{B[0, 1] n 7[0, oo) = 0, tq > 1}. 

The Beurling estimate gives us that the first summand is bounded by cr^^^. Thus, we only 
need to bound the second summand. To do this we let Ak be the event 

Ak = {Tk - Tk+i > 3"^'=+^) and tj - r,+i < 3"(^+^), for all j > k }. 

Note that 

r{rfc - Tfc+i > S-^'^+i)} = P°{ro - ri > (4/3)73} < Ciexp{-C2 (4/3)^}. (3) 

The last estimate is a standard estimate for Brownian motion and implies that F{Tk — Tk+i > 
3-(fc+i) i o. } = 0. Since Xlfcli 3"^^ < 1, we see that the event {tq > 1} is contained in Lik>oAk 
up to an event of probability zero. Using the Beurling estimate, we see that 

P^{fi[O,l]n7[O,oo) = 0,ro>l} 

oo 

< J]P^{fi[O,l]n7[O,oo) = 0,A,} 

oo 

< J2 ^'{B[0, Tk+i] n 7[0, oo) = 0, Tfc - Tfc+i > 3-('=+i)} 

k=0 

°° 1/2 

^ E<=(Fi^) ■fMr.-r.«>3-"«)}<cr'/=. 

k=0 

The last inequality uses 

□ 

Corollary 5.4. There is a c such that for every N, A and 9 < 2; there exists a coupling 
of the Brownian loop soup restricted to D [D] and a (l/N)-random walk soup restricted to 
3[D] such that the one-to-one correspondence of Theorem holds except on an event of 
probability at most c(A + 1) log A^^"^^/^)^ [respectively, c(A + 1) log A^^^/^) j^2-{b/i)ey 

Proof. It follows from Proposition 15.11 that the probability that a realization of the Brow- 
nian loop soup has at least one loop of time duration greater than cN^~'^ staying in ©, 
but not in for e = c(logA^/iV) is 0(A log A^^"'-^/^-*^). For general domains we get 
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0(A(log A^)*^^/^-* A^^"*-^/^^^) upon using Proposition 15.21 Therefore, if we consider the cou- 
phng of a (l/A^)-random walk soup and a Brownian soup and we restrict to those loops 
in a domain D, then we get the one-to-one correspondence of the loops as before, except 
on an event of probability c (A + i) iv2-(3/2)e logiV if D = D or an event of probability 
c (A + 1) iV^~(^/'*)^ logiV^/^ for general simply connected D contained in the unit disk. 

□ 

6 The dyadic approximation 
6.1 Introduction 

In this note we give a proof of the "dyadic" strong approximation for random walk bridges by 
Brownian bridges using the methods in Let Xi, X2, . . . be independent random variables 
with P{Xj = 1} = ¥{Xj = -1} = 1/2 and let S'n = Xi + ■ ■ ■ + X„. For positive integers n, 
let L„ = {2; G Z : P{S'„ = z} > 0}. If ^ G L„, {S'm'^''}^=o "^i^^ denote a process with the law 
of {Sm}m=o conditioned so that S'„ = z. 

We start with a definition. Suppose Z is a continuous random variable with strictly 
increasing distribution function F and G is the distribution function of a discrete random 
variable whose support is {ai,a2, . . .}. Then {Z,W) are quantile- coupled (with distribution 
functions [F, G)) if W is defined by 

W = aj if rj_ < Z < rj, 

where rj_, Vj are defined by 

F(r,_) = G(a,-), F(r,) = G(a,). 
The quantile-coupling has the following property. If 

F{ak -x) < G{ak-) < G{ak) < F{ak + x), 

then 

\Z -W\ = \Z - ak\ < X on the event {W = ak}. (4) 
We will need the following lemmas about the random walk; we will prove them in ^6.41 

Lemma 6.1. There exists eo > such that for every bi > there exist < c, a < 00 such 
that the following holds. Let N be a N{0, 1) random variable. For each integer n > 1, each 
integer m with \2m — n\ < 1, and every z E Ln, let 

Z = ^K",-) = !I^z+ L{1--)N, 
n \ n 

so that Z ~ N{^z,m{l - f )). Let W = be the random variable with C{W) = 

C{Sm'^^) that is quantile- coupled with Z. Then if \z\ < ton and ¥{W = w} > 0, 

E[e»l^-^l I = ^] < c exp (^^^^^^±^1 . (5) 
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Remark. For simple random walk, it is easy to show that (0) holds for eQii < \z\ < n (with 
perhaps different a,c), so it follows that the result holds for all \z\. However, we state the 
lemma only for \z\ < eon because this is all that we use. 

Lemma 6.2. There exist C2,&2;eo such that for every integer n > 2, every integer m with 
\2m — n\ < 1, every z & with \z\ < e^n, and every w E Tj, 

P{5„ = ^\S^ = z}< C2 n-'/' exp ! -b,^-^^^^-^^ 

I n 

Remark. We can actually show that this holds for any eo < 1 (with the constants C2, &2 
depending on eo). 

Let B denote a Brownian bridge, i.e., a Brownian motion in M conditioned so that 
Bq = Bi = Q (see ^6.21 for a more precise definition). If 2:1, 2:2 G M, n > 0, then 

r/-'-'-) := ^B„r. + "^,, + 1,,^ 0<t<n, (6) 

n n 

is the Brownian bridge conditioned so that Bq = zi, Bn = Z2. We write y^^"'^-* for y}"''^'^\ If 
{S^^'^\ B) are defined on the same probability space, we define 

A(n, z) = A(n, z, ^("-^^ B) = sup - '^^l. 

0<t<n 

(Recall that 5'^^"'^'' is defined for noninteger t by linear interpolation.) In Section we will 
prove the following. 

Theorem 6.3. For every b > 0, there exist < c,a,a < 00 such that for every positive 
integer n, there is a probability space on which are defined a Brownian bridge B and the 
family of processes {S^^'^^ : z G L„} such that 

E[e"^("'")] < cn"e^l"l'/". (7) 

Using Chebyshev's inequality, we get the following theorem as a corollary. 

Theorem 6.4. For every b > there exist < c, a < 00, such that for every postive integer 
n, there is a probability space on which are defined a Brownian bridge B and the family of 
processes {5"*^"'^^ : z G Ln} such that for all r > 0, 

F{A{n,z) > rc log n} < en""" e''"'/". 
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6.2 Brownian bridge 

If Wt denotes a standard one- dimensional Brownian motion, then the process 

Bt = Wt-tWi, 0<t<l, 

is called a Brownian bridge (conditioned so that Bq = 0, Bi = ). It can also be characterized 
as the continuous Gaussian process -Bt, < t < 1 with 

E[Bt] = Cov[fi„ Bt] = E[BsBt] = s{l-t), < s < t < 1. 

More generally, if Bt is a Brownian bridge and 

t — Si 



Xt = VS2- SiB,t-si)/(s2-si) +Xi + ( ) {X2 - Xi), 

S2 - Si 

is the Brownian bridge conditioned so that Xs^ = xi^Xg^ = X2- It is the continuous Gaussian 
process Xt, Si < t < S2 with 

E[Xt] = xi + (^^) ix2 - xi), Cov[X,Xi] = - - t) Si<S<t<S2. 

S2 — Si S2 — Si 

Lemma 6.5. Suppose B, B are independent Brownian bridges and N is an indepedent 
N{0, 1) random variable. Suppose < s < 1, and define Xt, < t < 1 by 

Xs = ^/sil -s)N 

Xt = ^rsBt/s +-Xs, 0<t<s, 
s 

. 1 -t 

Xf = VI - sB{t~s)/{i-s) + Xs, s<t<l. 

1 — s 

Then Xt is a Brownian bridge conditioned so that Xq = Xi = 0. 

Proof. This can be easily checked using the Gaussian characterization of Brownian bridges. 
The formulas are not mysterious. What we are doing is defining Xj, < t < 1, by first 
choosing X^ (using the appropriate distribution on X^), then defining Xj,0 < t < s, and 
-^t, s < t < 1, as appropriate Brownian bridges. □ 

We will need the following easy estimate for Brownian bridges. Let 

M = sup \Bt\. 

0<t<l 

Then there exist c, u such that for all a > 0, 

Ele"^^] < ge""'. (8) 

If Bt is replaced by a Brownian motion Wt, this estimate is standard using the reflection 
principle. That argument can easily be adapted to establish (jHl), perhaps with different c,u. 
(In fact, the maximum for Brownian motion stochastically dominates M so (jH)) holds with 
the same c,u, but we will not need this stronger fact.) 
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6.3 Proof of Theorem lOl 

It suffices to prove the result for b sufficiently small. We fix positive b < ^2/37 where 62 is 
the constant from Lemma lf).2l We let eo be the smaller of the two values of eo in Lemmas 
OandO 

In this proof, by an n-coupling we will mean a probability space on which are defined a 
Brownian bridge B and the family of processes {S"*-"'^-* : z & L„}. 

Note that for any n-coupling, if 2; G L„, St = S'f''^\ and Yt = y}^'^^ as in ©, then 



A(n, 2;)= sup \St — Yt\< sup + sup |Fj|<n + 

Q<t<n 0<t<n 0<t<n 



n + ^/n sup l-B- 

0<i<l 



Hence 



E[e"^("'")] < e2'^"E[exp{av/7^[sup \BtW}] < ce 



~ {2a+ua^)n 
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0<t<l 

where m, c are as in (jH)). Clearly, there exists oq = ao(6) > such that if a G (0,ao), then 
2a + ua? < 6 eg. 

Therefore, for any n-coupling inequality ((Tj) will hold with c = c, a = and a G (0, Oq) 
for all z & Ln with |z| > neo- For the remainder of this section, we will assume a < qq. We 
wil also assume that a is sufficiently small so that (0) holds with bi = 6/20. We now fix such 
a value of a, and we will show how to construct the n-couplings so that ^ holds for some 
c, a. 

We will use an induction. Clearly, we can choose n-couplings for n < 2 such that 

^[gaA(n,.)]g-fe|.|Vn Vz G L„, U < 2. 

We can assume without loss of generality that C > 1. We will show that there exists a 
constant c (which without loss of generality we can assume is greater than both 1 and the 
c in (jH))) such that for every positive integer s, if there exist n-couplings for all n < 2** such 
that 

^^aA{n,z)^^-b\z\yn ^ ^^^^^ ^g) 

then there exist n-couplings for all n < 2*+^ such that 

^[gaA(n,.)|g-fe|.p/n < cC{s). (10) 

The theorem follows easily from this claim. 

In order to prove the claim above, let 2** < n < 2*^^. We will show how to construct a 
probability space on which are defined a Brownian bridge and a family of processes {S'*^"'^-' : 
z G Ln, \z\ < neo} satisfying (fTUI) . Once this is done, we can adjoin, possibly after enlarging 
the probability space, the processes for \z\ > neo. Since c > c, firUj) will hold for these 
processes also. Hence, we assume \z\ < eon. For notational ease we will assume that n is 
even and we write n = 2k. Note that k is an integer with 2'^~^ < /c < 2*. (If n is odd we 
write n = k + {k + 1) and do a similar argument.) 

We define the n-coupling as follows: 
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Choose two independent fc-couplings 
satisfying 

Let N ~ iV(0, 1) and define the translated norma/ random variables = y^n/AN + ^. 
Define the quantile-coupled random variables as in Lemma lU^Tl Assume, as we may, 
that all these random variables are independent of the two fc-couplings chosen above. 
Note that a has been chosen sufficiently small so that (0) holds with bi = 6/20; i.e., 

E[e'^\^'-^'\ I W^ = w]<c exp[A !f!!±i!]. 

Let 

'''=UA~^^, + il-t)N \<t<l. ^''^ 
By Lemma inini Bt is a Brownian bridge. 



Let Sj^''^ = W, and 



Sl!''^'^ < m < fc, 

W + 5^^^r^'^ k<m<n. 



What we have done is to first choose the value of Sj^'^^ from the appropriate conditional 
distribution of 5*^ given Sn = z and then to define the other values of Sm'"^^ from the 
conditional distribution of Sm given Sn = z, Sk = . 

This gives us our coupling; we need to show that it satisfies (fTUj) . Note that 

\Z' - VT^I + max { A(A;, W\ S^^'''^^\B') , A{k, z - W, ^2(M-w^^)^ | ^ 
Therefore, 

Here we have used the the fact that our two fc-couplings satisfy (jHl) and the simple inequality, 

]g[gmax{Zi,Z2}j < +]E[e^2]_ Therefore, 

^^<cj:nW =-}exp|-^^|exp|6 1, 
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Since + < 5 maxji;;^, {z — w^} and k = n/2, this sum is bounded by 

9 max{w^, (2; — tLi)^} 



. ^ f 9 

c ^F{W' = w} expl-b- 



n 



We now split this sum into two pieces: |w — || < \z\/Q and jti' — || > \z\/Q. Iflw — 1| < \z\/Q, 
then max{w^, {z — w)^} < (2z/3)^; hence, 

P{H/^ = w;} exp j - & i — '-^ ^ \ < c exp{ }. 

To handle the other piece we make use of Lemma f6. 21 Since 62 > 376, we have 

iw - 

f>{W' = w]= ¥{Sk = w\ Sn = z}< cn-^'"^ exp{-376^ ^}, 

n 

and so, the sum over the second piece is bounded by 

c y n^l- exp{-376 ~ } exp{g " }. 

|«.-||>|^|/6 

For |w — || > |z|/6, we have [w — |)^ > ^ max{w^, (2; — w)^}. Hence, the sum is bounded 
by ^ 

c Vn-^/^ exp{-— — }, 
^-^ lo n 

which is clearly bounded by a constant. Therefore, 

That is, our n-coupling satisfies equation (fTUj). 
6.4 Local central limit theorem 

We will derive the local central limit theorem for conditioned simple random walks which 
is essentially a normal approximation for hypergeometric random variables. Although this 
is standard, we derive it here because the size of the error term is important for us. Our 
starting point is Stirling's formula with error 

n\ = v^n"+(^/2) g-n ^ 0{n-^)] . (12) 
Lemma 6.6. Suppose l,m,j are integers with m > 0, |/| < m/2, |j| < m/8. Then 
P{52™ = 2j + 21 I S^^ = 4/} = 



e^P <! ... + O(-) + 0(4) 



27r m(l — (//m)^) |^ 2m(l — (//m)^) m rrfi' 
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Proof. Throughout this proof we assume that |/| < m/2, \j\ < m/8. Without loss of gener- 
ahty we may assume that / > 0. Let 

/' 2m, W 2m \ 

p(2m, 2j I Al) = nS2m = 2j \ S,^ = Al} = ^^^^'-^^ 



f Am \ 
\2m+2l) 

Then, 

p(2m,2/ + 2(j + l) I 4/) =p(2m,2/ + 2j I 4/) ^ ' 



(m + j + 1)2-/2^ 
and if j > 0, 

(m - - /2 



p(2m, 21 + 2j I 4/) = p(2m, 21 \ Al) J] f \ J U2 n 

Aim + 2i — 2m — 1, 



i=0 

j 



p{2m,2l I Al) 



i=l 



(m + i)2 — /2 



By symmetry, p{2m, 21 - 2j \ Al) = p{2m, 21 + 2j \ Al). 
Stirhng's formula (fT^ gives 



f 2m\ 
\m+l) 



2 



p{2m, 21 I 4/) = = 2 



m 



{2I721) ^2nmil-{l/m)^) 
This gives the result for j = 0. To finish the proof of the result we need to show that 

Vl 4im + 2i-2m-l (2j)2 1 n/jl] 

^ (m + 2)2-/2 J- 2m(l-(//m)2)+^m^^^m3^' 



(13) 



Note that if / < m/2, j < m/8 and 1 < i < j, then 



Aim + 2i-2m-l 2 

< < -. 

- (m + i)^-P -3 



There exists Ci such that 



Therefore, 



log(l -X) +X+ yl < ClX^ 0<x<2/3. 



Aim + 2i — 2m — 1 , 



(m + i)^ — /2 

i=l ^ ' 



, j , V ^ Aim + 2i — 2m — 1 1 sr^Aim + 2i — 2m — 

(m + 2)2 - /2 (m + i)2-/2 ^ • 
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Note that 

1 1 2im 

(m + if - /2 " m2 - /2 " (^2 _ /2)2 + ^^>- 



Hence, 

Aim + 2i-~2m — l 1 4zm — 2m Sz^m' 



7 ^ = O— + ^ +> 5 



2 



Also, 



m + i)2 — /2 777, 7773 _ p Z-^ f^2 _ 12\2 

' i=\ i=l ^ ' 

^^1 ^ ^ 2j(j + l)m-2jm (8/3)jV 

m — (m^ — /^^ 

^,1 j\ 2j2m (8/3) j^m^ 

= 0( h — ) H ^ ' ' 

m rrv' ra^ — P [m? — Py 

1 J^Aim + 2i -2m-l,o j\ (8/3)jW 

- > ( ) = 0( h — ) + . 

2fr^ (m + i)2-/2 ^ (m2-/2)2 



The resuh follows immediately. □ 

One can deduce the following more general case from the result above. 

Lemma 6.7. There is ac such that for n >2 an integer, andm an integer with \2m—n\ < 1. 
// \z\, \w\ < c/n; z & Ln and w + G Lm- Then 

Tfl 

^{Srn = —Z + w\Sn = z} = 

n 



w"^ 1 w'^ 



2T:{m — [m? /n)){l — [z/nY) \ 2{m — {nP /n)){l — [z/nY) n rfi 

We now state without proof an easy large deviation estimate that follows from large 
deviations for binomial random variables. 

Lemma 6.8. There exists an r] > such that, for any a > 0, there exist C = C{a) < 00, 
and 7 = 7(a) > 0, such that for all z with \z\/n < rj 

Tfl 

^{\Sm z\ > amlSn = z}< Ce"^'" 

n 

Lemma 16.21 follows easily from Lemmas 16.71 and 16.81 



6.5 Coupling of conditioned distribution and normal 

In the remainder of this section we prove Lemma f6. 11 Note that we only need to prove the 
lemma for n sufficiently large. In order to simplify the notation we will assume that n is 
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even and hence m = n/2. If is odd, one can do the same argument. We will use a slightly 
weaker form of Lemma 16.71 more precisely, we will assume that 



exp 



w 



'n 



\w\ 



(14) 



where cr^^ = (n/4) [1 — (z/n)^]. We have replaced the 0{l/n),0{w'^/n^) terms with the 
larger 0{l/y^),0{\w\^/n'^) terms, respectively. Our reason for doing this is that more 
general random walks satisfy local central limit theorems with this weaker error term, and 
it is useful to know that the arguments in this section only use the weaker form. 
Let X denote a A^(0, 1) random variable, let 



Wn^z be the random variable with distribution of S^^'^^ that is quantile-coupled 



and let W 

with X. Note that W is also quantile-coupled with Z and Z. Let us write F = Fn^z for the 
distribution function of Z and G = Gn.z for the distribution function of W. It follows from 
Lemma 16.91 below that there exist c, e and such that for all n > N, for all z E Ln with 

\z\/n < e, and all w with \w — {z/2)\/n < e, 

{w - {z/2)f 



F w - c 1 



n 



-)) < G{w-) < G{w) < F{w + c(l 



{w - iz/2)f 



n 



))• (15) 



It follows from Q and p5|) that the quantile-coupling satisfies 

{W-{z/2))' 



W <c 



1 + 



n 



for all n > N, provided that \z\, \W — {z/2)\ < en. Also 



1 



n 



z 



< c- 



n 



provided that 1^1, \W — {z/2)\ < en (which implies the estimate \Z\ < d n). Therefore, for n 
large enough the quantile-coupling satisfies 



W\<c 



^ ^ {W-{z/2)f ^ z^- 
n n 



for|z|, \W - {z/2)\ < en. 



For any e' > 0, straightforward exponential estimates show that there exists an a > such 
that (0) holds for \w\ > e' n. Hence to get Lemma lOl it suffices to prove the following 
estimate. We have written the estimate for random variables without the odd/even parity 
issues of simple random walk and hence have dropped a factor of 2. If we have a random 
variable that satisfies (Hlj), i.e., that is supported only on even or only on odd integers, we 
can convert it to a random variable on all the integers by dividing the mass at k equally 
between k and k — 1. 
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Lemma 6.9. For every c, e there exist ci,ei,Ni such that the following holds for every 
positive integer n > Ni and every G [1/c, c]. Suppose S is an integer random variable 
such that for every integer \j\ < i n, 

nS = j} = -j^== exp + , (16) 

where 

m\<c — + ^ 



Assume additionally that for every a there exist positive a, b such that f{\S\>an}< ce"''", 
where c and b do not depend on S . 

Then if G denotes the distribution function of S and F denotes the distribution function 
of an N{0, a'^n) random variable, 



F\^-Ci[l + —]j<G{x-l)<G{x + l)<F\^ + Ci[l + —]j, \x\<ein. (17) 

Proof. This is a straightforward estimate of sums and integrals. It suffices to prove (fTTjl for 
integer x, and by symmetry we can assume x > 0. Since we only need to establish the result 
for n large when we write inequalities in this proof we will only be asserting that they are 
valid for n large enough. 

Let F = l-F, (5 = 1-G', <l = 1-$, where $ is the distribution function of a A^(0, 1). 
Let q{y) = qiy;n,a^) = {27ra^n)-^/^ e-y'/^^'''''\ 

We will deal with the somewhat easier case < \/3dy/n, at once. Note that (fTHjl implies 
c'/ y/n < G{x) — G{x — 1) < c'/ y/n for |x| < \f^yfn. By definition, 

G{X) = J2nSn =J} = Y1 ^(3) + J2i^{Sn = j} - q{j)]. 

j>x j>x j>x 

Also, 



5^g(j)=0(^)+ / q{y)dy. 



Straightforward estimates using the assumptions give 



J2\nsn=j}-q{j)\ = o{^). 

■yii' 

j>x 



Hence, we see that (fT7|) holds for \x\ < \p^^n. 

It remains to prove ()17p for x > v^V^- To make our strategy more intuitive for the 
reader we state now the simple fact about the standard normal distribution that lies at the 
heart of our lemma. The interested reader is refered to [51 for many more details, and further 
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reading on the KMT approximation. For any A > 0, there exist a c and an e > such that 
for cr^/n < z < en, 



e^^$(^)<$(^-c^), (1^ 
awn oJn on'^i'^ 



e"^^$( — -=) > $( — -= + c 



Assume now that y^v^ < 3; < n^^^. The choice of 5/8 as the exponent is somewhat 
arbitrary. One could take any exponent strictly less than 2/3 as the argument below shows. 
We have 

= j} - q{j)\ < 0(exp{-C3nV3}) + c ^ q{j). 

j>x j>x 

Using that x is an integer and that y^q{y) is decreasing for y > ^/Sdy/n, we see that 

foo 3 1 1 poo 



But for s > 1, 



3 e-^'/2 



oo poo 

Therefore, 

2 2 S 

j>x 

Also, using simple estimates we obtain 



Hence, 



^-^ oJn n 



X , r ^,X^ 



G(x) = <l>(^)exp{0(-)}. 
a^n n'^ 

The result for x in this range now follows from (fTSj) . 

Assume now that n^^^ < x. Note that there is a constant c, depending only on c, such 

that q{y)e~^ is decreasing for x < ?/ < 2cn, and such that 

q{y)e~^dy < / q{y)e~^dy. 



Ix-l Jx-2 

From this we see that there is an e2 such that if a; < £2^, then 



G'(x) < / q{y)e~^ dy + G{cn - I) < I q{y)e~^dy< I q{y)e~^dy. 



2x 

x-1 Jx-2 
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Hence, 

G(x) < — ^ . 

a^/n 

Similar arguments can be used to obtain 

G{x) > $( — 1=). 

The result for x in this range follows from ()18|1 . and this concludes the proof of the lemma. 

□ 
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